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Error  Probability  Characteristics  for  Multiple  Alternative 
Communication  With  Diversity,  but  Without  Fading 

Introduction 

The  M-ary  character  error  probability,  Pc,  for  orthogonal  multiple  alternative 
communication  with  D-fold  diversity,  was  evaluated  in  reference  1  for  a  wide  range 
of  parameter  values.  It  was  presumed  there  that  the  slow  signal  fading  on  the  D 
diversity  channels  was  independent  and  that  the  received  signal  energy  was 
distributed  exponentially  (Rayleigh  fading  of  received  signal  voltage). 

!'  is  of  interest  here  to  reconsider  the  system  performance  for  the  case  where, 
although  the  receiver  was  designed  for  D-fold  diversity,  there  is  in  fact  no  signal 
fading.  The  fractionalization  of  the  received  signal  into  the  diversity  channels  then 
results  in  poorer  performance  than  had  the  signal  been  confined  to  one  channel  and 
processed  coherently.  This  situation  can  arise  naturally  in  practice,  as  for  example, 
as  a  result  of  multipath  arrivals,  or  it  ean  come  about  intentionally  in  the  system 
design.  It  can  also  be  of  interest  in  testing  a  multiple  alternative  hardware  design 
under  controlled  laboratory  conditions,  where  the  time  or  cost  of  simulating  actual 
fading  conditions  is  excessive. 

The  basic  framework  and  background  for  the  communications  technique  con- 
idered  here  have  already  been  presented  in  reference  I  and  will  not  be  repeated,  for 
t  ie  sake  of  brevuy.  The  reader  is  presumed  to  be  familiar  with  the  time-bandwidth 
duration  and  separation  constraints  listed  in  the  above  reference,  particularly  pages 
2-4  and  appendix  B. 
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Definitions  and  Technical  Results 


The  source  selects  one  of  M  equi-probable  symbols  and  transmits  a  signal  on  D 
diversity  channels  (time  and/or  frequency)  to  a  receiver.  The  D  channels  devoted  to 
each  of  the  M  alternatives  are  disjoint  (in  time  and/or  frequency)  with  each  other 
and  with  the  channels  for  the  other  alternatives  (see  reference  I).  The  receiver 
employs  matched  filtering  on  each  of  the  totality  of  MD  cells  of  interest.  For  each 
signal  alternative,  the  corresponding  D  envelope-squared  matched  filter  outputs  are 
sampled  appropriately  in  time  and  summed.  The  largest  of  these  M  decision 
variables  is  then  declared  to  be  the  transmitted  signal  alternative. 

The  additive  noise  at  the  receiver  is  assumed  to  be  white  Gaussian  over  the  total 
band  of  the  possible  received  signals,  with  a  (single-sided)  power  density  level  of  N(, 
watts/H/.  It  is  shown  in  appendix  A  that  the  M-ary  character  error  probability  of 
this  communications  technique  is  given  by 


P  =  1 
e 


(I) 


Here  “deflection”  statistic  d,  is  given  by 


d  £ 

_T  __  _T  total  received  signal  energy 
2  Nq  noise  power  density  level 


(2) 


where  ET  is  the  total  received  signal  energy  on  the  D  channels;  ln(x)  is  the  modified 
Bessel  function  of  order  n  and  argument  x;  auxiliary  function  Bn(x)  is  defined  as 


=  f  *  «,s 


'0 

and  (reference  2,  equation  6.5.1 1 ) 


,00  =  Z  rr  *k 

1  k=0  K- 


(3) 


(4) 


is  the  exponential  power  series  through  the  term  x". 

It  is  interesting  and  worthwhile  to  observe  from  (I)  and  (2)  that  the  exact  frac- 
tionali/ation  of  the  total  received  signal  energy  E,  amongst  the  D  diversity  branches 
is  immaterial  in  so  far  as  system  performance  is  concerned.  The  available  signal 
energy  E ,  could  be  divided  equally  in  the  D  branches,  or  it  could  be  concentrated  in 
just  one  branch;  in  either  ease,  error  probability  Pc  is  identical  (assuming  M  and  I) 
are  unchanged). 

l  or  M  =  2,  expression  (I)  can  be  evaluated  in  closed  form.  Several  alternative 
representations  for  this  binary  error  probability  are  given  in  appendix  B.  along  with 
an  efficient  program  for  this  special  case.  The  program  for  evaluation  of  (I)  in 
general  is  given  in  appendix  C. 
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(Graphical  Results 

For  a  particular  value  of  M,  the  character  error  probability  Pt.  in  (I )  is  a  function 
of  the  number  of  diversity  channels,  D,  and  the  total  energy-density  ratio  (signal-lo- 
noise  ratio)  E,/Nil.  We  have  plotted  Pc  on  a  normal  probability  ordinate  versus 
Er/N(1  on  a  dB  abscissa  (i.e.,  10  log  ET/N((),  with  D  as  a  parameter,  in  the  range  10  15 
to  10  1  for  P,.  Figures  1-1 1  correspond,  respectively,  to  M  =  2,  3,  4,  8,  16,  32,  64, 
128,  256,  512,  1024.  The  curve  for  D  =  1  in  each  figure  corresponds  to  the  NO 
FADING  results  in  reference  I,  figures  1-7. 

We  have  discovered  no  simple  rule-of-thumb  for  the  additional  signal  energy 
required  to  maintain  a  fixed  P  as  D  is  increased,  which  covers  the  whole  range  of  I) 
and  \1.  Nor  did  we  find  a  simple  expression  for  Pt.  w  hen  it  is  small  10  '),  despite 
considerable  effort.  The  best  we  can  do  is  to  observe  that 

Pg  as  (M  -  1)  Pg2  for  small  Pg  ,  (5) 

and  then  use  result  ( B- 1 0)  for  the  binary  error  probability,  Pc„  D  is  arbitrary  in  (5). 
Further  development  of  this  approach  appears  in  appendix  B. 
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Et/N0  in  dB 


Figure  11.  M-ary  Character  Error  Probability  for  M  =  1024 
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Correct  Decision  with  a  Threshold 


The  processing  system  above  presumed  that  one  of  the  M  signal  alternatives  was 
always  transmitted.  A  more  general  situation  occurs  when  the  additional  ease  of  no 
signal  present  at  all  is  allowed.  Then  the  largest  of  the  M  decision  variables  is 
compared  with  a  threshold.  There  are  three  types  of  probabilities  of  interest  for 
signal  present  (reference  3,  pages  2-5);  they  are  the  probability  of  a  missed  decision, 
the  probability  of  an  incorrect  decision,  and  the  probability  of  a  correct  decision. 
The  Iasi  is  given  by  a  slight  generalization  of  (1)  to  (reference  3,  equation  5) 


where  V  is  the  threshold  value. 


f  or  signal  absent,  we  then  have  false  alarm  probability  (reference  3.  equation  2) 


fa  -  1  ’  Probjmax(ylt  ...  ,  yMJ  <  V  for  dy  =  oj 


D-l 


=  1 


-  Jo  dyy  mp  (~y2/^ 


M 


=  1  -  [bo.jIv2^)]  .  (7) 

Here  we  used  the  limit  of  the  integrand  of  (6)  as  dj-'O  and  (3).  Result  (7)  agrees 
with  reference  3,  equation  25,  as  it  must,  since  both  receivers  are  processing 
identical  noise-only  channel  outputs.  No  numerical  results  on  this  more  general 
situation,  (6)  and  (7),  are  presented  here. 

Extension  to  Fading  Signal 


All  the  earlier  results  in  this  report  have  pertained  to  constant  amplitude  signal 
components  on  the  D  diversity  branches.  Now  we  presume  that  there  is  slow  signal 
fading,  meaning  that  the  total  received  signal  energy  F(  is  a  random  variable.  In 
particular,  we  let  energy-density  ratio  (signal-to-noise  ratio) 


d  .  h. 

Kt  -  2 


have  probability  density  function 

.v 


3(rt)" 


Rt  exp^-l^/p) 

uv+1  r(v  ♦  1) 


for 


"T  > 


(S) 


<l» 


w*--  1 
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»T  VNo 

V  +  1  ~  V  +  1 


,  V  >  -  1  . 


Here  ET  is  (he  average  received  total  signal  energy  on  the  p  branches,  and  v  is  a 
constant,  not  necessarily  an  integer. 

This  very  general  model  of  fading  was  presented  and  utilized  in  reference  4,  pages 
1 1-12,  in  the  investigation  of  a  maximum  likelihood  detector.  For  example,  the  four 
signal  fading  cases  considered  by  Swerling  (reference  5)  are  subsumed  by  (9)  for 
particular  choices  of  v;  see  table  I . 

Table  I.  Particular  Cases  of  General  Fading 


V 

0  D-l  1 

2IJ-I 

Swelling  Case 

1  2  3 

4 

Cases  I  and  3  correspond  to  common  signal  fading  on  a ||  D  brandies  whereas 
cases  2  and  4  correspond  to  independent  and  identically  distributed  signal  fading  on 
each  of  the  D  branches.  In  particular,  case  2,  v  =  D  -  I,  corresponds  to  ex¬ 
ponential  probability  density  function 

Pl‘R)  =Tr"p("'5T)forR>0  ’  •  (III 

for  the  received  energy-density  ratio,  E,/N0,  on  one  branch;  this  js  Rayleigh  fading 
of  the  received  signal  voltage  amplitude  on  each  branch.  Equation  (9)  is  a  D-fold 
convolution  of  (II)  when  v  =  D  -  I.  Case  4,  v  =  2D  -  I,  corresponds  to 
probability  density  function 


Pl«)  -=!««*("£<>) 


r.5 


for  R  >  0  ,  , 


for  the  received  energy-density  ratio,  E(/N0,  on  one  branch.  Equation  (9)  is  a  D- 
lold  convolution  of  (12)  when  v  =  2D  -  1. 

The  average  character  error  probability  is  evaluated  in  appendix  D;  it  is  given  by 

-  1  -  [(D  -  1)!  (I>  *  1)V*1J 

•  f0  dt  ,D'‘  iFi(v  * 1;  rh-  *)  •  (in 

This  single  integral  can  be  evaluated  numerically  for  any  v  of  interest;  of  course,  M, 
D,  and  E,/Ntl  need  to  be  specified  also  for  this  numerical  procedure. 

As  a  special  instance  of  fading  (9),  consider  case  2  in  more  detail: 
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V  =  D  -  1, 


Rt  Et/Nq 

P  =  -g-  =  — p —  =  average  energy-density  ratio  per  branch.  (14) 
Then,  if  we  use  reference  2,  equation  13.6. 12,  and  (3)  above,  ( 1 3)  specializes  io 


P  =  1  - 
e 


/*“  r  rx  im_i 

io  dt  q'(t)  [Jo  *  q»<X,J 


(I'l 


where 


.D- 


Ct)  = 


1  f„,  t » o  . 


(D  -  1)!  (y  ♦  l)1 


D-l  ,  , 

x  exp(-x)  -  ^  - 

qo(x)  =  '  (5  -l) :  for  x  >  0  '  (16) 

are  the  probability  density  functions  of  the  decision  variables  formed  by  the  sum  of 
H  envelope-squared  filter  outputs  for  signal-present  and  noise-onlv.  respectively. 
Result  (15)  js  identical  to  refcrcnu  1,  equation  C-l.  It  is  what  we  would  have  gotten 
had  we  averaged  the  signal  probability  density  function  in  (I)  —  the  function 
multiplying  the  bracket  —  with  respect  to  the  signal  strength,  before  expressing  P 
in  integral  form.  In  fact,  a  more  general  integral  than  (15)  has  already  been  en¬ 
countered  and  evaluated  in  reference  3,  equation  26,  where  a  thresholding  operation 
was  included  in  the  receiver  processor  for  the  fading  signal. 

l  or  M  2,  general  fading  result  (13)  can  be  evaluated  in  closed  form;  these 
special  cases  are  presented  in  appendix  D,  especially  (D-5)  and  (D-8). 

Discussion  ahd  Summary 

t  his  report  has  addressed  the  problem  of  numerically  assessing  the  additional 
signal-to-noise  ratio  (i.e. ,  energy-density  ratio)  required  when  a  nonfading  signal  is 
broken  into  D  (unequal)  components  and  combined  incoherently.  Some  related 
stuJies  into  the  cost  of  imperfections  of  receivers,  or  the  cost  of  lack  of  knowledge 
of  the  received  signal  structure,  are  presented  and  evaluated  in  references  6  and  7. 

Since  the  performance  of  the  processor  considered  here  for  the  nonfading  signal 
depends  only  on  the  total  received  signal  energy,  regardless  of  how  Iractionali/ed  it 
may  be  among  the  diversity  channels,  these  results  also  apply  to  a  situation  where 
uncertainly  in  the  time  of  arrival  or  the  doppler  shift  of  the  received  signal  causes 
the  receiver  to  search  over  several  channels,  even  though  the  particular  received 
signal  alternative  occupies  only  one  channel.  The  D-fold  summation  for  each  M  at  v 
alternative  leads  to  a  loss  in  performance;  stated  alternatively,  the  cost  of  the  un¬ 
certainty  is  additional  received  signal  energy  required  in  order  to  maintain  the  same 
quality  of  performance,  P  . 

I  or  a  large  signal-to-noise  ratio,  the  M-ary  character  errot  probability  is  given 
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approximately  by  (5)  in  terms  of  the  binary  error  probability,  P.2.  The  usual  ap- 
proaeh  is  then  to  use  (A-3)  with  M  =  2,  or  to  use  (B-4),  to  evaluate  Pc,.  However, 
the  situation  frequently  arises  where  the  characteristic  functions  of  the  decision 
variables  can  be  evaluated  fairly  easily,  whereas  the  corresponding  probability 
density  functions  and/or  cumulative  distribution  functions  cannot.  In  that  ease,  an 
alternative  representation  of  Pc2  directly  in  terms  of  the  pertinent  characteristic 
functions  would  be  useful.  This  problem  is  addressed  in  appendix  E.  with  the 
following  results  for  the  binary  error  probability: 

Pe2  "I ffo<E)fl(-« 

■  -rarjT  ¥  • 

+ 

’7*?/0  f  •  (17, 

Here  fDU)  and  f,(4)  are  the  characteristic  functions  of  the  noise-only  and  signal- 
plus-noise  decision  variables,  respectively.  C_  is  a  contour  in  the  complex  |-plane 
along  the  real  axis,  with  a  small  downward  identation  at  |  =  0;  C+  is  a  similar 
contour  indented  upward  at  4  =  0.  The  last  form  in  (17)  can  be  particularly  useful 
for  numerical  evaluation  of  Pe2  for  some  characteristic  functions. 

Another  alternative  to  evaluation  of  M-ary  error  probabilities  is  to  use  bounds 
which  are  simpler  to  compute  than  the  exact  result.  This  is  attractive  if  the  bounds 
are  tight,  at  least  for  the  range  of  signal-to-noise  ratios  and  error  probabilities  of 
interest.  Some  general  bounds  on  error  probability  (which  were  not  used  here)  are 
presented  in  appendix  F  lot  a  situation  that  includes  interference  as  well  as  noise. 
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Appendix  A 

Derivation  of  M-ary  Character  Error  Probability 


To  avoid  duplication  o I  effort,  we  will  rely  heavily  on  reference  8,  pages  2,  3,  8,  9 
and  appendixes  A  and  B,  especially  (B-7MB-12).  We  obtain  directly  the  probability 
density  function  of  the  signal  decision  variable  (sum  of  D  filter  outputs)  as* 


and  that  for  any  one  of  the  M-l  noise-only  decision  variables  as 

po(y)  =  I  Sd"-3!)';  exp  (-  2)  for  y  >  0  •  (A  2) 


A  correct  decision  is  yielded  only  if  x  is  greater  than  all  M-l  independent  noise 
variables;  the  probability  of  this  event  is 


P 


c 


dx  px(x) 


dy 


po(y) 


M-l 


( A  -  3 ) 


But  the  integral  on  y  in  (A-3)  and  (A-2)  is  given  by  (3)  as  B0  ,(\/2).  Then  by  use  of 
(A- 1 )  and  the  substitution  x  --  y2,  (A-3)  becomes 


P 

c 


exp 


(- 


4/2) 


,D-1 


f 


dy  y 


-y2/2 


-i(dTy)[BD-i(y2/2)] 


M-l 


(A -4) 


I  lie  M-ary  character  error  probability  is  P.-i-  Pc. 

Several  checks  on  (A-4)  are  available.  First,  for  the  trivial  case  M  =  I,  (A-4) 
yields  P  =  1,  which  agrees  with  the  fact  that  there  are  then  no  noise-only  channels 
to  worry  about  at  all.  Second,  as  d  ,“*0,  (A-4)  reduces  to  -•  I /M;  this  agrees  w  ith 
the  observation  that  we  have  a  completely  random  selection  with  zero  sigual-10- 
noise  ratio.  Third,  for  D  I ,  (A-4)  agrees  with  reference  1 ,  equation  9A  (when  the 
latter  is  corrected  for  the  typographical  omission  of  the  factor  cxp(-F. ,  /Nt,)). 

A  program  for  the  evaluation  of  P  .  =  1  -  Pc  via  (A-4)  is  given  in  appendix 


*A  “deflection”  interpretation  of  d^  is  available  from  reference  8.  equations  (B-10),  (B-7), 
(B-4),  (A-4),  and  (B-5). 


A  I  A-2 
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Appendix  B 

Binary  Error  Probability 


This  appendix  will  deal  exclusively  with  i he  special  ease  of  M  2,  a  binarv 
decision.  Then  (I )  becomes,  with  (he  identifications  in  (3)  and  (4), 


To  evaluate  the  integral,  we  used  reference  9,  equation  6.631  I;  the  final  trans¬ 
formation  to  (B-l )  employed  reference  2,  equation  13.1.27.  An  alternative  ex¬ 
pression  to  (B- 1 )  is  available  via  reference  2,  equation  13.6.9: 


P 


e2 


1 


D'l 
k=0  2* 


CD-I) 

k 


(B-2) 


If  we  expand  the  polynomial  in  d^/4  that  occurs  in  (B-l )  and  (B-2),  and  interchange 
sums,  we  obtain 

exp(-d^/4)  D-l  (dbi)  D-l-n  (D  +  n)  . 

p  =  — LI  L  y  LJ_/  y  —■ - 1  . 

6  T  n=0  2n  n!  jTo  2J  j !  (B-3) 

All  three  expressions  above  lor  P.-,  are  finite  sums  of  positive  quantities  and  are 
reasonable  for  computer  evaluation,  even  for  fairly  large  1). 

Another  expression  for  P  -,  -  1  -  P  2  is  available  b\  recalling  ( A- 1 )  and  (A  2) 
and  changing  (A-3)  to 

Pc2-/0  d".W  Jy  d»PIW  ■  „,41 

But  from  (A-l),  upon  making  the  substitution  \  t:  and  using  lelerenee  10. 
equation  I ,  we  find  for  the  inner  integral  in  (B-4), 
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D-l 


.2  .2' 

t  *  dT 


/•  /•  /  t  + 

jy  dXP>W  ’/i/2dt  ''r\-—2 

•  «(dT.  r1'2)  ■ 


>-i(V) 


(B-5) 


Diet)  (B-4)  yields 


c2 


/; 


D-l 


dy  i  _  saLxm  q  (a  y>/2) 

2  (D  -  1)!  UV  T  ' 


L 


=  [20'1  (D  -  1)  -•]  du  u2D_1  exp(-u2/2)  QD (dT ,  u)  .  (B.6) 

Ai  this  point,  we  will  use  (he  new  integral  result 
dx  x2N  1  exp(-p2x2/2)  Q  (a,  bx) 


2N_1(N  -  1)! 
2N 


M+N-l 


1 


-A 


N-l/M+N-l' 


2 

D  +  b 


’  Sol 


eN-l-n  (A) 


w  here 


A  =  " 


2  2.2 

p  +  b 


(B-7) 


(B-8) 


I  bis  integral  and  several  other  new  integrals  of  QM  functions,  along  with  their 
derivations,  will  be  presented  fairly  soon  by  this  author  in  another  NUSC  technical 
i  eport . 

I  he  use  of  (B-7)  and  (B-8)  in  (B-6)  yields 
exp 


e2 


I  a  (*m 


(B-9) 


I  xpanding  the  partial  exponential  expansion  in  powers  of  d-j/4  according  It'  (4), 
and  interchanging  summations,  we  obtain  for  the  binary  error  probability 


e2 


exp(-d2/4)  D-l  (d2/4)  D-d-k  (2 D-l 

2D  - 1  ~  jjTt  2-»  \  n 

2  k=0  * 1  n=0  '  n 


(B-IO) 


I  his  result  is  very  similar  to  (B-3),  and  of  course  yields  identical  numerical  values; 
however,  a  direct  verification  of  the  transformation  that  would  lake  the  inner  sum 
m  (B-IO)  to  the  inner  sum  in  (B-3),  or  vice  versa,  has  not  been  discovered. 
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(B-10)  is  a  very  efficient  form  for  numerical  evaluation,  because  the  summand  ot 
the  inner  sum  (on  n)  does  not  depend  on  the  index  (k )  of  the  outer  sum;  this  ad¬ 
vantageous  property  is  not  true  of  (B-3).  Thm  for  a  given  I),  wc  can  compute  the 
partial  sums  on  n  in  (B-10)  by  starting  at  k  =  I)  -  1  and  merely  adding  one  ad¬ 
ditional  term  to  the  inner  sum  for  each  downward  unit  step  in  k;  these  D  partial 
sums  are  all  stored  and  then  used  for  the  outer  sum  starting  at  k  -  0.  Thus,  what 
appears  to  be  a  double  sum  in  (B-10)  actually  can  be  computed  by  means  of  two 
single  sums  of  D  positive  terms.  Furthermore,  the  inner  sunt  on  n  in  (B-10)  can  be 
done  once  for  a  value  of  D  of  interest,  and  then  used  repeatedly  for  different  values 
of  d-j.  A  program  incorporating  these  observations  is  listed  at  the  end  of  ibis  ap¬ 
pendix.  All  of  the  results  thus  far  apply  to  arbitrary  values  of  energy -density  ratio 
E,/N0. 

One  reason  for  dwelling  on  the  evaluation  of  the  binary  error  probability  is  that 
for  large  signal-to-noise  ratio,  i.e.,  large  energy-density  ratio  F|/Nit.  the  M-ary 
character  error  probability  is  given  approximately  by 

P  »  (M  -  1)  P  _  for  small  P  .  (B  ID 

e  e2  e 

I  bis  result  may  be  seen  to  be  valid  when  it  is  realized  that  the  probability  ot  any 
particular  noise-alone  decision  variable  exceeding  the  signal  decision  variable  is  a 
rare  event  when  P  is  small.  Thus  the  probability  of  two  or  more  noise  variables 
exceeding  the  signal  variable  is  very  rare  indeed,  and  wc  need  only  account  lot  the 
single  noise  variable  case.  Since  these  events  are  disjoint  and  there  are  M-l  such 
variables,  (B-J I )  follows  immediately. 

For  large  d}/2  =  Er/N0,  an  asymptotic  form  for  Pc2  can  be  developed  from 
(B-10);  express 

exp(-d^/4)  (dJ/4) 
e2  "  22D-1  (D  -  1)  ! 


1  +  2D (D  -  1)  (D  -  1)  (D  -  2)  (2D  -  D  ♦  1) 

2  ,  „  /  _  \  2 


d‘/4 


Therefore,  asymptotically. 


(4/4y 


(B  12) 


exp  (“^7/4)  (dx/4) 


D-l 


Pe2  22D-1  (D  -  1) !  3S  dT  ^  “  ' 


(B-1 3) 


But  reference  to  (B-12)  reveals  that  in  order  for  (B- 13)  to  be  a  good  approximation 
to  P  ,,  rather  than  just  an  asymptotic  result,  the  requirement  d,  »  l)  must  be 
satisfied.  This  requirement,  however,  is  typically  too  large  a  signal-to-noise  ratio 
requirement  to  be  practically  useful.  Wc  have  found  it  necessary  to  resort  to  (B-IO) 
for  the  binary  error  probability  for  the  range  of  values  of  interest  here. 
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A  program  for  the  binary  error  probability  Pi2  in  BASIC'  for  the  Hewlett- 
Packard  9845  is  presented  below.  It  utilizes  (B-10)  et  seq.,  as  discussed  above.  The 
variables  are  self-explanatory  and  consistent  with  the  earlier  text. 


I.  1  Fir.i'o  Pc  tor  D-fold  diversity,  but  no  fadirvj.  TP  6472 
1C  D»12  !  D>0  NUMBER  OP  DIVERSITY  CHANNELS 

20  Et  1  E»  n0>0  TOTAL  SIGNAL  ENERGY  - S I NGLE-3 1  BED  NOISE  DENSITY 

SO  D 1 =  D -  1 

40  B  2  =  B  *  2 

50  A=Et  ri0  2 

60  DIM  S  Y 1 00 z 

70  REDIM  SsDl > 

SO  S  Y  D  l  >  =  T  =  1 

SO  FOR  J  =  1  TO  D 1 
100  T  =  T*  Y  D2- J  >  J 

no  s  *  d  i  - 1  > = s  <:  d  -  j  > + t 

120  NEXT  T 

ISO  s*S*'0> 

140  T=  1 

1 60  FOR  K = 1  TO  D1 
1  r. 0  T-  T* A  - K 

170  S=itT*S<K> 

ISO  NEXT  K 

ISO  Re  =EXR  Y  -  A  '•  *S -'2  B2-l> 

20o  PRINT  D  *"  5  0,  "Et  ri0  ="  ;  Et  r,0,  "P«  »•*  ;Pe 
210  END 


=  20 


F' «  =  3.  4 5131 9 3 3  1 3 ?£-0 3 


E »  r,  0 


[  R  6473 


Appendix  C 

Program  for  M-ary  Character  Error  Probability 


The  expression  for  Pc  is  given  by  (1).  Since  we  cannot  integrate  to  we  terminate 
the  integral  at  y  =  L.  The  error  incurred  in  doing  this  is  given  by 


exp 


dy  yD  e"y  /2  ID-i(dTy)  [BD-i(y2^)] 


< 


dy  y 


«M°~2  exp[-  i(L  -  di)2] 

\V  (2tt)1/2  (L  -  dT)  ' 

We  used:  the  fact  that  Bn(x)  is  upper-bounded  by  I,  as  may  be  seen  from  (3);  the 
asymptotic  behavior  of  ln(x)  for  large  x  as  given  by  reference  2,  equation  9,7.1;  and 
an  integration  by  parts  to  yield  the  final  form.  For  given  values  of  D  and  d,  (M  has 
dropped  out  of  this  error  bound),  we  step  !  up  by  units  of  1  from  the  value  d,  +  I 
until  (C -1 )  is  sufficiently  small. 

We  chose  error  10  12  in  the  enclosed  program.  The  reason  that  the  integral  in  { I  > 
must  be  evaluated  very  accurately  is  that  we  must  subtract  the  integral  result  front  i 
in  order  to  gel  the  error  probability.  If  we  want  to  evaluate  Pc  accurately  in  the 
range  of  10- \  then  we  need  the  integral  to  8  or  9  decimal  places  to  counteract  the 
effect  of  subtraction  from  I  required  by  form  (1). 
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We  numerically  integrate  (1)  for  y  in  the  interval  (0,  L)  via  the  Trapezoidal  rule, 
with  automatic  halving  of  the  interval,  until  a  stable  result  occurs.  It  is  shown  in 
reference  1 1,  pages  55  and  58,  that  the  Trapezoidal  rule  is  excellent  for  integrands 
w  hose  odd  derivatives  at  the  end  points  of  integration  are  equal.  For  the  application 
here,  the  integrand  of  (I)  may  be  shown  to  behave  as  y2l)M  as  y-*-0  +  .  Thus  a  high 
number  of  odd  derivatives  are  zero  at  y  =  0  +  and  virtually  zero  at  y  =  L  for  large 
DM,  meaning  that  the  correction  terms  to  the  Trapezoidal  rule  vanish.  In  fact,  (I) 
was  also  integrated  via  Simpson’s  rule  and  via  Simpson’s  rule  with  end  correction 
(reference  12,  pages  414-418);  neither  performed  as  well  as  the  Trapezoidal  rule, 
w  hich  is  incorporated  in  the  program  for  Ptf  below.  Because  of  the  very  small  values 
for  the  exponentials  and  the  very  large  values  for  the  Bessel  function,  it  was  found 
necessary  to  first  evaluate  the  logarithm  of  the  integrand  of  (1)  and  then  ex¬ 
ponentiate  it  prior  to  its  use  in  the  Trapezoidal  rule. 


1  ' 

M - ar  y  Re 

t  or 

D -tol 

d  d  i  " 

€  r  £  1  t.  y  , 

,  b 

ut 

no  f  4-^; 

i  r't  9  .  T  K 

£473 

10 

11  =  S' 

i 

ft  >  1 

NUMBER 

OF 

SI 

GNAL  ALTERNATIV 

fc'S 

20 

D  =  1  2 

i 

DO 

NUMBER 

OF 

III 

VERS  IT 

CHANNEL 

S 

?  0 

Et  <10  =  20 

i 

E»  r.O 

0 

TOTAL  ( 

SIGNAL 

ENERG' 

'SINGLE- 

SIDED 

4  0 

Dt =S0T ■ 

2  *  E  t 

n0  > 

i 

TOTAL  1 

‘DEFLE 

CT  I  ON 11 

STATIST  I 

c 

50 

PRINT  " 

M  =" 

;  M,  "  D 

=  -;d 

,  " Et. nO 

=  " 

;  Et 

nO , "Dt 

=  " ;  Dt 

SC' 

COM  D , D 

1  ,  Ml 

,Dt  ,E 

t  nO 

?e> 

M 1  =  M  -  1 

so 

D 1  =  D  -  1 

so 

3 1  =0 

1  00 

S  2  =  D  t 

lie 

A  =  D  -  .  5 

120 

B= . 4  D« 

ft 

1  so 

32  =  -1 2+1 

140 

C  =  S  2  -  D  t 

150 

E  =  B*S2'V 

ft  *  E  X 

P  f  -  .  5 

*  C  *  C  ) 

zC 

1  SO 

IF  E  > 1 E 

-12 

THEN 

130 

1  70 

—  K  'z?  “ 

S 1  )  * 

•  5 

ISO 

34  =  2 

1  SO 

55  =  s FNS 

s  S 1  > 

♦  FNS  < 

S2  >  >* 

.5 

200 

S  £  =  0 

210 

FOR  S7  = 

1  TO 

S4-1 

STEP 

2 

iCO  S6«S6*FNS<S1«-S3*S?) 

2-30  NEXT  3? 

2-0  PRINT  US  INC  "M.  12DE,7D";  1 -(.S5+S6  >*S3,  S4 
250  S 3=33*. 5 
2cC'  54  =  S  4  *2 

270  35=S5+S6 

2*0  GOTO  2O0 
2*0  1 


NOISE 


DENSITY 


t-2 


! 


it*  it*  if)  u~'  u“'  u-'  it*  i»*'  u“> 
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36i©  DEF  FNS  •  V  1  ‘  INTEGRAND 

310  COM  D, D1 , Ml , Dt , Et n0 
330  T  = . 5  *  V  * V 

330  e=i-e::P‘  -t  wfnect.di  > 

340  IF  B '  =  0  THEN  RETURN  0 

350  S=B*L0G‘ Y  Dr  >-T+FNI  nxl  og<  Dt  *Y,  D1  )+M  1  *i_G0< B  > -Et  riO 
■36 0  RETURN  EXP<S>*Dt 

370  FHEND 

330  i 

330  DEF  F NE  s  X ,  N >  !  SUM  OF  X'K-k!  FROM  K  =  0  TO  N 

400  3  =  T  =  1 

410  FOR  k  =  1  TO  N 

430  T  =  T  *  X  K 

430  S  =  3  +  T 

440  NEXT  1 

450  RETURN  S 

460  FNEND 

470  ! 

4 3 0  D E F  F N I  t‘i x  1  c* 9 f-  X «  H  )  !  LOG  OF  I  rt  C  x 

430  IF  -0  THEN  520 

500  IF  N=0  THEN  RETURN  1 

10  IF  N >0  THEN  RETURN  0 

20  A=.5*X 

30  F  =  1 

40  FOR  1=2  TO  N 
50  F  =  F*I 

60  NEXT  I 

70  F=N*LOGCfl 3-LOGCF3 
30  A=A*Fi 

90  S=T=1E-S? 

600  FOR  1=1  TO  500 
6  1  61  T  =  T  *  A  /  <  I  *  s  N  +  I  )  > 

620  S=S  +  T 

630  IF  T * 2 E  1  1  =5  THEN  660 
640  NEXT  I 

650  OUTPUT  O;"5O0  TERMS  AT  “;X;N 
6  6  O  I  n  x  1  •:>  g = F + L  0  G  <  £ .  >  +  2  O  0 .  3  2  4  9  0  3  0  9 

670  RETURN  Inxlog 
630  FNEND 


M  -  3  D  =  12 

- . 22 3 3 30444766E+O 1  2 

620  1602  -1 19 30E  +  00  4 

. 423327566O0OE-02  8 

.  *846 3 ; 5 1 0  900E - 0 1  16 

. 1S42167O25O0E-01  32 

. 1 342 1 66 1 O50OE-0 1  64 

. 18421661 170OE-O1  128 


Et nO  =  20 


C-VC-4 
Reverse  Blank 
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Appendix  I) 

Average  Error  Probability  for  Fading  Signal 

The  energy-density  ratio  variable  RT  was  defined  in  (8),  and  its  probability 
density  was  given  in  (9)  for  a  fading  signal.  From  (I),  the  error  probability  expressed 
in  terms  of  R  r  is  given  by 

•*p(-Rt) 

D-l 

(2rt)  2 

(D-l) 


The  average  error  probability  is 


Pe  '  f’  dRT  P(*t)  %(Rt)  ' 

Substitution  of  (9)  and  (D-l )  in  (D-2)  and  interehange  of  integrals  lead  to 


(D-2) 


P  =  1 
e 


-  c « s  [h-iMT1 


•'O  x-y-  y  r(v  +  1) 

W 


(D-3) 


In  the  inner  integral,  let  RT  =  1 2/2  and  use  reference  9,  equation  6.631  I.  Sim¬ 
plification  leads  to 

1-1 

-1 , 


P  =  1 
e 


T-xD-l^  ,x,  ,  ,,v+ll  f°°  .  2D-1  -y2/2 

-  |^2  (D  -  1)!  (y  +  1)  J  J  dy  y  e  J  / 


(D-4) 


Use  of  the  variable  t  =  y2/2  in  (D-4)  yields  the  result  already  quoted  in  (13).  Results 
(D-4)  attd  (13)  hold  for  general  v  and  M.  When  v  is  chosen  as  the  special  value  D-l, 
the  resulting  average  error  probability  is  given  by  ( 15)  and  (16). 

Specialization  to  M  =  2 


Now  we  let  parameter  v  in  (9)  be  general,  but  we  set  M  =  2.  Thus  we  will  develop 
the  average  binary  error  probability  for  a  general  signal  fading  model.  Upon  use  o' 


D-l 
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(3)  and  (4),  (13)  yields  for  M  =  2, 

-.-1 


Pe2  .  [(0  -  1) 


!  (U  •  lH 


g  &  /;  *  ‘D-**k  ,',(v  *  o 

*  1)D'*'V  V  f(d  -  1  -  v  D  +  k;  D;  — — 

(F  +  2)°  ft,  '  k  A»  *  V  \  •  "  * 


D-l-v  D-l 


(U  +  2) 


This  finite  sum  of  Gaussian  hypergcometric  functions  is  not  too  useful  in  general. 
However,  for  the  special  case  v  =  D  -  1  considered  in  (14)-(16),  (D-5)  immediately 
yields 

—  1  v-1  /D-l+k\/u  +  l\k  , 

62  (u  +  2)°  to  \  k  AU  +  2>  (D-ft) 


this  last  result  can  also  be  derived  directly  from  (15)  and  (16)  for  M  -  2,  by 
reference  to  (3)  and  (4). 

A  much  more  appealing  result  than  (D-5)  for  average  binary  error  probability  P  , 
for  general  v  is  attained  if  we  start  with  the  binary  no-fading  result  in  (B-9)  and  use 
identification  (8): 


„  exp  (-Rt/2)  D-l  /2d_A  D-Un  , 

-  '  220'1  nMn/  i~g  k'  \  2  /  •  (I 

Then  the  average  binary  error  probability  is,  by  use  of  (D-7)  and  (9), 

Pe2  =  fQ  dR T  P(Rt)  Pe2 

.[a20-1  ,«*  rev  .  I)]’'  E(2Dn-‘)  °t"  -4-f 

L  J  n=0v  n  '  k=0  k!  2K 

■  /o  dRp  ^  "p  (■  T  •  ■ 

r,2D-i  i.  U\U*1T‘  Dfr  fo.i\  D-frn  (v  *  1Jk  /  m  \ 

T  r J  SI-)  £  -n—  [ttw) 


U  "  2 
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This  closed  form  result  is  a  finite  sum  of  positive  quantities  and  holds  for  general 
values  of  t/(>—  1 ).  The  summations  are  very  similar  to  those  encountered  in  non¬ 
fading  binary  error  probability  (B-10).  A  program  for  the  evaluation  of  (0-8), 
which  takes  advantage  of  the  observations  made  under  (B-IO),  is  given  below.  I  or 
v  =  0-1,  (D-8)  must  reduce  to  (D-6);  we  have  not  discovered  the  transformation 
that  accomplishes  this  corroboration,  but  we  have  confirmed  the  equality 
numerically  for  v  =  D  -  1,  as  well  as  for  the  alternative,  more  general  result  (0-5). 
No  numerical  results  for  (0-8)  have  been  presented  here. 


1  1  Average  binary  P«  for  D-fold  diversity  with  fading. 

10  DM2  !  D>0  HUMBER  OF  DIVERSITY  CHANNELS 

20  Mu*6.1  <  No >  - 1  MEASURE  OF  FADING,  EOS.  9-10 

JO  Mo* 3.7  1  Mu>*0  MEASURE  OF  SIGNAL-T0-N0I3E  RATIO 

40  PRINT  "D  *";D,”Nu  =" J  Nu, "Mo  *“;Mu 

50  D 1  =  I'  -  1 

60  D  2 - D  *  2 

70  A-Mu-  '<2  +  Mu) 

SO  DIM  S< 100 > 

90  REDIM  SCDl > 

100  S  >.  D 1'  * T  =  1 

HO  FOR  JM  TO  D 1 

120  T  =  T *  ■ .  D 2 - J  >  r  J 

130  SsDl- J>=S<D- J>+T 

140  NEXT  J 

1  SO  S  =  S  < O ) 

160  T  =  1 

170  FOR  K* 1  TO  D 1 

180  T  =  T  *  < N u  +  K  1  *  A  r K 

190  S  =  S  +  T*SOO 

200  NEXT  K 

2  1  0  Pe  =  S  ‘  2'  D2  - 1  >  ■*  <  1  + .  5*Mu  >  A  <  Nu+ 1  >  > 

220  PRINT  "F'e  =";Pe 

2  JO  END 


b  =  1.2  Mu  *  6.1  Mu  =  3.7 

F  t  =  6  .  1 0  C 1 5  5 1 6  8 6  3  E  -  O  3 


D-3/D-4 
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Appendix  K 

Binary  Error  Probability  Via  Characteristic  Functions 


l.et  p„(y)  and  p((x)  be  the  general  probability  density  functions  of  a  noise-only 
decision  variable  and  a  signal-plus-noise  decision  variable,  respectively.  The  binary 
probability  of  error  is  then 


Pg2  =  Prob(x  <  y) 


dx  p:(x) 


dy  pQ(y) 


(ID 


Now  by  reference  13,  equation  6,  the  exceedance  probability  in  ( (• - 1 )  is  given  in 
terms  of  the  characteristic  function  f  ,(4)  according  to 


dy  p0(y)  =  1  -  p0M 


1  1 
2  *  i2ir 


f  X  fo(^  exp(-i^x) 
*4oo  ^ 


=  15?  Jc  ~f  V0  exPC'iW  •  (E-2) 

The  integral  with  a  slash  on  it  is  a  principal  value  integral;  the  contour  C_  is  a 
contour  on  the  real  axis  of  the  complex  {-plane  except  for  a  small  downward 
identation  at  {  =  0.  Employment  of  (E-2)  in  (E-l)  yields 

/+oo  - 


■  12?  /  f  fo®  fl<-« 


-fk  L  ffo(-«  £.«>  •  (HI 

where  C  ,  is  indented  upward  at  4=0.  Birth  contours  C  ,  start  at  4  =  and  go 

ll>4  -  +  CO, 

Expression  (E-3)  gives  binary  error  probability  Pe2  directly  in  terms  of  charac¬ 
teristic  functions  fQ({)  and  f,(|),  by  means  of  a  contour  integral.  Since  characteristic 
functions  satisfy  the  relation 

f(-£)  =  f*(0  for  real  £,  ,  (1-4) 


L-l 
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an  alternative  real  integral  form  for  (E-3)  is 


_  4.  CO 

Pe2=H  1  f  Im{£o(0  f: 


E-2 


(-O) 
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Appendix  F 

Bounds  on  Error  Probability 

We  consider  the  situation  where  there  are  three  different  variables  upon  which  a 
decision  as  to  the  transmitted  multiple  signal  alternative  must  be  made.  1  lies  are. 
respectively,  the  signal,  noise,  and  interference  variables.  In  particular,  lei  random 
decision  variable  S  denote  the  output  of  the  signal  channel;  S  may  itself  he  the  sum 
of  several  diversity  branch  outputs.  Let  N  be  the  output  of  the  maximum  of  the 
noise  channels  outputs,  of  which  there  may  be  many;  each  noise  channel  output 
may  itself  be  the  sum  of  several  diversity  branch  outputs.  And  let  I  be  the  output  of 
the  maximum  of  the  interference-plus-noise  channel  outputs,  of  which  there  may  he 
many;  each  interference  channel  output  may  itself  be  the  sum  of  several  diversity 
branch  outputs. 

For  independent  signal,  noise,  and  interference  channel  outputs,  the  probability 
of  a  col  l  ect  decision  is 

Pc  =  Probfl  <  S,  N  <  S)  =  /dx  ps(x)  P.(x)  Pn(x)  ,  (I  I) 

where  p,(x)  is  the  probability  density  function  of  random  variable  a,  and  I’  ,(\)  is  its 
corresponding  cumulative  distribution  function.  We  also  express 

P.(x)  =  Prob(I  <  x)  =  1  -  Q. (x)  , 

Pn(x)  =  Prob(N  <  x)  =  1  -  (^(x)  ,  0  -2) 

where  Q  ,(x)  is  the  exceedance  probability  of  random  variable  a. 

The  probability  of  character  error  is 

Pe  =  1  '  Pc  =/dx  Ps00!1  ’  {1  '  Vx)Hl  ’  Qn(x)}] 

=  fdx  Ps(x)[q.(x)  +  (^(x)  -  Qi(x)  Qn(x)] 

=  Prob(I  >  S)  +  Prob(N  >  S)  -  Prob(I  >  S,  N  >  S)  . 

(1-3) 

The  last  quantity  in  (F-3)  is  generally  quite  small  for  practical  cases  of  interest.  But 
in  any  event,  we  have  the  upper  bound 

Pg  <  Probfl  >  S)  +  Prob(N  >  S)  (and  Pg  <  1)  .  (14) 


In  addition,  since 
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